FERMIL AB-Pub-94 /OO 1- A 
astro-ph/9405036 
submitted to Physical Review D 



MASSIVE TAU NEUTRINO 
AND SN 1987A 

G. Sigl 1 - 2 and Michael S. Turner 1 - 2 ' 3 

1 Department of Astronomy & Astrophysics 
Enrico Fermi Institute, The University of Chicago, Chicago, IL 60637-1433 

2 NASA/Fermilab Astrophysics Center 
Fermi National Accelerator Laboratory, Batavia, IL 60510-0500 

3 Department of Physics, Enrico Fermi Institute 
The University of Chicago, Chicago, IL 60631-1^33 

ABSTRACT 

The emission of MeV-mass tau neutrinos from newly formed neutron stars is consid- 
ered in a simple, but accurate, model based upon the diffusion approximation. The tau- 
neutrinosphere temperature is found to increase with mass so that emission of massive tau 
neutrinos is not suppressed by the Boltzmann factor previously used, (m^/T^) 1,5 exp(— m u /T u ). 
where T u ~ 4 MeV — 8 MeV. If the tau neutrino decays to electron neutrinos, then for short 
lifetimes [t v < 10~ 3 sec) the location of both the tau and electron neutrinospheres can be 
affected, and, for very short lifetimes (t v < 10~ 6 sec) its temperature falls below 1 MeV, 
in conflict with neutrino observations of Supernova 1987A (SN 87A). Using our results, we 
revise limits to the mass/lifetime of an MeV-mass tau neutrino based upon SN 87A. Our 
constraints, together with bounds based upon primordial nucleosynthesis and the laboratory 
mass limit of around 30 MeV, exclude the possibility of a tau neutrino more massive than 
0.4 MeV if the dominant decay mode is radiative. Finally, we speculate on the possible role 
a 15 MeV — 25 MeV tau neutrino might play in the supernova explosion itself. 



1 Introduction 



In most extensions of the standard model neutrinos are massive and unstable. A massive 
neutrino can decay through standard electroweak interactions into "visible" channels such 
as v — > z/7 or v — > z/77, and if it is more massive than about 1 MeV, v —>■ v e e + e~ and 
v — > v e e + e~^j. It can also decay through new interactions into "invisible" channels, e.g., 
v — > v'cf) or f — ► 1/ + z/V, where is a massless (or very light) boson (e.g., majoron or 
familon) . 

Many of the most stringent limits to the radiative decay of a massive neutrino come from 
the nondetection of a gamma-ray burst coincident with the neutrino burst [|TJ observed from 
Supernova 1987A (SN 87A) |@» 0» B @' S 01 an< ^ f rom the fact that the energy observed in 
supernovae in the form of bulk kinetic energy of the expanding envelope and light (about 
10 51 erg in total) are far less than that carried away by a neutrino species (about one third of 
the the binding energy of a neutron star or 10 53 erg) ||. These limits are based on neutrino 
fluxes calculated within the framework of standard supernova theory where it is assumed 
that all three standard model neutrino species are massless and stable. The upper limits 
to the masses of the electron and muon neutrinos, about 7 eV and 250 keV respectively, are 
much less than the temperature of their neutrinospheres, on the order of 4 MeV and 8 MeV 
respectively, and so this is certainly a good approximation. The current laboratory limit 
to the tau-neutrino mass is slightly more than 30 MeV || and so such an approximation 
is not necessarily valid for the tau neutrino. Previously, the mass dependence of the tau- 
neutrino energy flux has been obtained by simply multiplying the energy flux of a massless 
neutrino species by a Boltzmann factor, implicitly assuming that the temperature of the tau 
neutrinosphere does not depend its mass. 



Because of the importance of the SN 87A constraints to the tau-neutrino ||10|| , we felt it 
worth considering more carefully the effect of finite mass and lifetime on the energy fluence 
of tau neutrinos. To this end, in Sec. 2 we develop a simple model of the transport of 
MeV-mass tau neutrinos based upon the diffusion approximation. In Sec. 3 we use the 
tau-neutrino energy fluence calculated in this model to re-derive the constraints that follow 
from observations of SN 87A. 



2 A Simple Model for Tau-neutrino Transport 

2.1 Basic considerations 

The energy transport within, and the cooling of, a newly born neutron star is predominantly 



due to neutrinos [|Tl|]. Because the neutrino mean free path is very short in a neutron star, 
it takes neutrinos a few seconds to diffuse from the center of the star to the surface. In 
considering the transport of tau and mu neutrinos the most important processes are elastic, 
neutral-current scattering off nucleons, vN —>■ uN, and pair creation/annihilation processes 
involving nucleons and pairs, <-> vv and NN <-> NN + vv. For electron neutrinos 
charged- current processes involving nucleons and electrons/positrons are also important. 
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We define the neutrinosphere to be the surface beyond which neutrino number is ef- 
fectively conserved, i.e., pair creation/annihilation processes have become ineffective and 
chemical equilibrium ceases to be maintained. Because the elastic-scattering cross section 
is significantly larger than the corresponding pair creation/annihilation cross sections mu 
and tau neutrinos continue to scatter off nucleons well beyond the neutrinosphere. However, 
elastic scatterings conserve neutrino number and approximately conserve neutrino energy, 
and so they have little effect on the neutrino energy flux. Thus, the flux of mu and tau 
neutrinos is determined by the conditions at the neutrinosphere. 

In passing, we note that the situation for electron neutrinos and antineutrinos is different 
since they can scatter inelastically through charged- current interactions, e.g., v e + e ± — > 
v e + e^, u e + n — > e~ + p, or v e + p — > e + + n. The electron neutrinosphere is very close 
to the mu/tau neutrinosphere, and in the massless limit the numbers of electron, mu, and 
tau neutrinos radiated are essentially equal. However, charged- current, inelastic scatterings 
that take place beyond the neutrinosphere serve to lower the kinetic temperature of electron 
neutrinos, bringing it closer to the local temperature, ultimately to around 4MeV. This 
explains an apparent paradox: namely, how electron neutrinos can be characterized by a 
temperature that is about a factor of two smaller than mu/tau neutrinos and yet have a 
fluence that is about the same, rather than a factor of (T Ue /T u ) 3 ~ 1/8 smaller. 

We shall use a simple model based the diffusion approximation to calculate the tau- 
neutrino energy flux; as mentioned above, diffusion is a very good approximation at the 
tau neutrinosphere. The location of the tau neutrinosphere involves an interplay between 
scatterings, pair creation/annihilations, and the local physical conditions (temperature and 
density). We assume that the underlying physical conditions in the nascent neutron star 
are the same as those in the standard case (three massless neutrino species). Since, energy 
transport by tau neutrinos only accounts for about one third of the total energy transport 
(and doesn't change radically for the masses of interest) this assumption is well motivated. 
Finally, in actually solving for the flux of tau neutrinos from the neutrinosphere we use a 
technique akin to the freeze-out approximation used in relic-abundance calculations in the 
early Universe; by comparison to direct numerical integration of the transport equations we 
find it to be a very good, time saving approximation. 

In the massless limit (m v < T^) we find that the temperature of the tau neutrinosphere 



is 7.5 MeV, in agreement with standard calculations Further, our results are very 

insensitive to the input parameters (cross sections, model for physical conditions, and so 
on), indicating that our various approximations have little effect on the energy flux of tau 
neutrinos we calculate. Both these facts give us confidence that our model provides reliable 
results. 



2.2 Interactions and physical conditions 

In the diffusion approximation the flux of massive tau neutrinos = —DVn where n is the 
neutrino number density and D = Xv/3 is the diffusion coefficient, given by the product 
of the neutrino mean free path A and velocity v divided by three. Throughout, all energy- 
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dependent quantities are understood to be averaged over a thermal distribution. Inside the 
neutrinosphere the temperature is given by the local temperature and outside the neutri- 
nosphere it is taken to be constant and equal to the temperature of the neutrinosphere. 

In general, we approximate thermally averaged quantities by a functional form that 
extrapolates correctly to the ultrarelativistic and nonrelativistic limits; for example, 



3T 



m v + 3T 



1/2 



Our results are sufficiently robust that the fact that quantities do not have the correct 
behaviour in the semirelativistic regime is of no consequence. 
The diffusion coefficient is given by 

D(r) = , (2) 

3n B ^ N 

where tib is the local baryon number density and ovn is the cross section for elastic scattering 
of massive tau neutrinos on nucleons (in a hot, nascent neutron star baryons exist in the 
form of free nucleons rather than nuclei; see Ref. | II). We approximate the elastic-scattering 
cross section by |12| 
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(3) 



depending on whether the massive tau neutrino is of the Dirac or of the Majorana type. 

We parameterize the physical conditions in the nascent neutron star during its early, very 
hot phase (first 5 sec — 10 sec) by 



n B (r) 



T(r) 




(4) 



where spherical symmetry has been assumed and for the parameters r , n , T and k we 
choose 

r = 30km, n = 5.97 x 10 35 cm" 3 , T = 5MeV, k = 5; (5) 

see e.g., see Ref. [ 1 1]. Our results are insensitive to the precise parameters chosen. 

Deep inside a hot neutron star pair creation/annihilation processes involving nucleon 
neutrino-pair bremsstrahlung dominates those involving pairs because of its strong density 
dependence and the fact that electrons are very degenerate. However, in the important region 
for our calculation, near the neutrinosphere, the two processes are of comparable importance. 
We write the total cross section for massive neutrino pair annihilation as their sum, 
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where the numerical constant a^N ~ 2.93 x 10 3 . The cross section for the process is 



from Ref. |y| and neglects the fact that electrons are mildly degenerate; the cross section 



for the nucleon bremsstrahlung process is from Refs. [13, 14] and is computed in the dilute- 
medium approximation which is justified for the highly subnuclear densities around the 
neutrinosphere. 

We also allow for tau-neutrino decays. The rate of neutrino decay and inverse decay is 
controlled by the neutrino lifetime and energy, with the rate for both processes is reduced 
by the Lorentz 7 factor in the usual way. We take 22.5 MeV to be the typical energy for tau 
neutrinos and approximate 7 by 



22.5 MeV 
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(7) 



In deriving the transport equation we assume that the daughter products of tau-neutrino 
decays are in thermal equilibrium so that detailed balance can be used to supply the inverse- 
decay rate. This assumption is certainly justified if the decay products are u e 's, u^s, 
pairs, or photons and applies to many of the decay processes of interest. 

Finally, for the equilibrium number density of a massive tau neutrino we take 



n, 
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(8) 



which extrapolates to the correct form in the ultrarelativistic and nonrelativistic limits. 



2.3 Neutrino transport equation 

The neutrino transport equation follows from adding pair annihilation/creation and de- 
cay/inverse decay to the continuity equation for the neutrino number density 

n + V0 = -o vD v{n 2 - n 2 J - {n - Weq) , (9) 

where <fi is the neutrino flux and we have used detailed balance to relate pair creation to pair 
annihilation and inverse decay to decay. 

The use of the diffusion approximation and the assumption of spherical symmetry and 
stationarity (n = 0) allows the transport equation to be written as an ordinary second-order, 
nonlinear differential equation 

Dn" + (d' + -d) n' = a vV v[n 2 - nl) + , (10) 

V r ) 4 777 

where prime denotes derivative with respect to r. 

To solve Eq. ( |TUD two boundary conditions must be specified. To arrive at the first we note 
that the diffusion approximation breaks down around the surface of last scattering (r = r^is) 
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where neutrinos begin to stream freely out of the star. For r ^> r^s the flux is therefore 
given by <fi « nv/2 where the factor of two comes from averaging over the angle between 
the propagation direction and the surface normal. On the other hand, where the diffusion 
approximation applies (r <C r^ff) we have <fi = —DVu. Matching these two relations leads 
to the first boundary condition 

-DVn| rdW =(n V /2)| rdB . (11) 

Now for the second; well inside the neutrinosphere the number density of tau neutrinos 
is very close to its equilibrium value because annihilation/creation interactions are occurring 
rapidly. We can therefore linearize Eq. (pf) by setting n = n eq in terms that are of zeroth 
order in the small quantity (n — n cq ) and solve for n to first order in the deviation from 
n eq . By so doing and solving for n(r) well inside the neutrinosphere we obtain our second 
boundary condition and Eq. fllQ|) can be integrated outward from this inner boundary. (By 
varying the position of the inner boundary we have shown that our results do not depend 
significantly on it, so long as it lies considerably inside the neutrinosphere.) 



2.4 The tau neutrinosphere 

While we have solved the neutrino transport equation for the flux of massive tau neutrinos 
by direct numerical integration, we have found that an approximation based upon the "freeze 
out" of the number of tau neutrinos provides a good approximation that requires far less 
computational time. In particular, the tau-neutrino energy fluence calculated agrees within 
about 50%. Given that we are using an approximate model for the nascent neutron star, 
this accuracy seems adequate. 

Recall, at the neutrinosphere the actual neutrino number density begins to differ from 
its equilibrium value because the mean free path for neutrino pair annihilation/creation 
becomes too large to maintain equilibrium.^ To find the neutrinosphere we algebraically 
solve Eq. (|10"D for the radius r ns where n — n cq becomes of order unity by setting n — n eq = n eq 
and n 2 
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From this, we approximate the power radiated in massive tau-neutrinos, 
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(13) 



1 The analogy with the freeze out of a massive particle species in the early Universe is a good one, and the 
approximate technique we use to solve the transport equation is very similar to the freeze-out approximation 
used in relic-particle calculations; see Ref. fl5|] . 

2 As it turns out, there is a small mass range where this procedure does not work well because the term 
on the left-hand side becomes zero. To insure that r ns is a smooth function of m v and t v we actually take 
the maximum of this term and 0.1 times the modulus of an analogous term with a minus sign between the 
second and the first order derivatives which could otherwise cancel each other. 
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where T v = T(r ns ) is the temperature at the tau neutrinosphere. While the power agrees 
well with the standard results in the massless limit, we shall actually use the ratio of the 
power for a massive tau neutrino to that of a massless species and the standard results for 
a massless neutrino species to obtain absolute results for the energy fluence (in any case, 
the total energy radiated depends upon the cooling time of the hot neutron star). The 
number fluence of and energy carried by a massless neutrino species (both neutrinos and 
antineutrinos) are about 1.4 x 10 10 cm -2 and 10 53 erg respectively. 

Our results for the tau neutrinosphere temperature and tau-neutrino energy fluence are 
shown in Figs. 1 and 2 as a function of tau neutrino mass and lifetime. In the limit of 
long lifetime (t > 1CT 2 sec) the tau neutrinosphere temperature remains constant at about 
7.5 Me V until a mass of about lOMeV, after which it rises steadily, reaching around 14MeV 
for a mass of 100 MeV. The reason for the rise is simple: As the mass increases beyond 
10 MeV the tau-neutrino equilibrium abundance at the "massless neutrinosphere" decreases 
exponentially; pair-annihilation processes vary as the neutrino abundance squared and can- 
not reduce the neutrino number density to its small equilibrium value for large tau-neutrino 
masses; thus, the tau neutrinosphere moves inward, to a higher temperature where the abun- 
dance is larger. The net result is that the tau-neutrino energy fluence relative to a massless 
neutrino species decreases far more slowly than it would if one (incorrectly) assumes that 
the neutrinosphere temperature is independent of neutrino mass; see Fig. 2. 

In Figs. 1 and 2 the effect of the tau-neutrino lifetime upon the neutrinosphere is also 
shown. For lifetimes greater than about 10~ 2 sec there is little effect; as the lifetime is 
decreased, which corresponds to increasing the rate of decays and inverse decays, the tau- 
neutrinosphere temperature decreases. The reason is simple: Inverse decays and decays 
become very effective at maintaining chemical equilibrium, further and further out. In 
fact, they are so effective that for tau-neutrino lifetimes shorter than about 10" 3 sec the tau 
neutrinosphere temperature drops below 7 MeV, and presumably, the electron neutrinosphere 
temperature is lowered too. For very short lifetimes the effects of decays and inverse decays 
are so severe that the background model is likely to be affected too, calling into question its 
validity. 

Finally, a brief, but important, comment on the robustness of our results. Consider the 
tau-neutrinosphere as defined above. Its location depends upon the model of the ambient 
conditions as well as the rates for elastic scattering and pair-creation/annihilation. Vary- 
ing the model for the ambient conditions, e.g., k and T , affects the tau-neutrinosphere 
temperature by only a few tenths of an MeV. Likewise, it depends very little upon the 
overall normalization of the scattering and annihilation rates adopted: In particular, it is 
the product of the rates that determines T„, and changing that product by a factor of ten 
only changes the tau-neutrinosphere temperature by about 1 MeV and the energy flux by 
less than a factor of four for long-lived (i.e. r v > 0.1 sec) tau neutrinos. Insensitivity to input 
parameters is also apparent in the small differences in the results for Dirac and Majorana 
neutrinos, where the product of the rates differs by a factor 0.1(m v /T u ) 2 for m v > 20 MeV. 
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2.5 Wrong-helicity neutrinos 

A massive Dirac neutrino has two additional helicity states, ur and z/£, which in the absence 
of new interactions only interact by virtue of the mismatch between chirality and helicity. 
For every ordinary process, e.g., i>l + N — > ul + N, there is a spin-flip analog, here, ul + N — > 
Vr + N, which is suppressed by a factor of {m u /2E v ) 2 . 

For masses much less than 1 MeV the suppression factor is very large and the mean-free 
path of a wrong-helicity neutrino is much larger than the radius of a neutron star. Wrong- 
helicity neutrinos that are produced simply stream out, tending to accelerate the cooling of a 
hot neutron star, e.g., the one associated with SN 87A, and thereby shortening the duration 
of the predicted neutrino signal. The duration of the neutrino burst observed by KII and 
1MB has been used to exclude a Dirac neutrino of mass between about 10 keV and 1 MeV 

For masses much greater than 1 MeV the suppression factor is not significant, and 
wrong-helicity neutrinos become trapped and are radiated from a "wrong-helicity" neu- 
trinosphere. Using our diffusion model we have estimated the location of the wrong-helicity 
neutrinosphere. (The dominant interactions for both scattering and creation/destruction are 
spin-flip neutrino-nucleon scatterings; see Ref. ||13|| .) For a mass of 1 MeV the wrong-helicity 
neutrinosphere temperature is about 40 MeV, and wrong-helicity neutrinos will quickly rob 
the core of the hot neutron star of most of its thermal reserves. As before, this is excluded 
by the KII and 1MB data. As the mass increases the temperature of the wrong-helicity 
neutrinosphere decreases, reaching a temperature comparable to that of the proper-helicity 
neutrinosphere for m v = 5 MeV. For m u ^> 5 MeV both neutrinospheres coincide, thereby 
roughly doubling the energy flux (compared to a Majorana neutrino). 

The two neutrinospheres are coupled (through spin-flip interactions); to properly take 
account of both helicity states of an MeV Dirac neutrino is beyond the scope of this study. 
Since we are only interested in neutrino masses greater than a few MeV, where wrong-helicity 
neutrinos carry off only about as much energy as proper-helicity neutrinos, to be conservative 
we simply ignore the energy carried off by wrong-helicity neutrinos. (Doubling our energy 
fluxes would change our results very little.) 

3 Constraints Based Upon SN 87A 

In this Section we re-examine the important mass-lifetime limits to an MeV-mass tau neu- 
trino based on SN 87A using the more accurate tau-neutrino fluence computed here. In 
general, our limits are significantly more stringent as the tau-neutrino energy fluence is 
larger than had previously been assumed. The limits we discuss are based upon three ob- 
servations: (i) the total "visible" energy, the kinetic energy of the expanding shell of matter 
(about 10 51 erg) and the optical light output (about 10 49 erg), was less than 10 51 erg [jnj; (ii) 
the gamma-ray fluence around the time when the neutrinos were detected was very small 
compared to the enormous neutrino fluence, F v ~ 10 10 cm" 2 Q; and (iii) 19 neutrino events 
were detected by the KII and 1MB detectors |I[ , consistent with antielectron neutrinos emit- 
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ted with a kinetic temperature of about 4 MeV |TI| . In deriving our mass-lifetime limits for a 
tau neutrino of mass greater than a few MeV, we follow closely the treatments in Refs. [j7|, |18[ 
and also assume that the decay rate is dominated by visible channels (i.e., daughter products 
include photons or pairs). 

3.1 Decays inside the progenitor 



The progenitor of SN 87A had a size of about 3 x 10 12 cm |L7]. For tau- neutrino lifetimes of 
the order 100 sec or smaller an appreciable fraction of the tau neutrinos emitted decay inside 
the progenitor and deposit energy. Since both pairs and photons have small mean free 
paths inside the progenitor and thus should deposit all their energy within the envelope of 
the progenitor, our arguments do not depend whether the daughter products are pairs or 
photons. The fraction of tau neutrinos that decay inside the progenitor 

/inside = 1 - exp(-t inside /r L ), (14) 

where inside ~ lOOsec/f is the time it takes a tau neutrino to escape the progenitor and 
t~l = JT U is the lifetime in the progenitor rest frame. 

We demand that the energy deposited inside the progenitor inside be less than about 
10 51 erg to be consistent with the observations of SN 87A (and other type II supernovae): 

inside « /insidelO 53 erg pj ) S ^ ^ (^) 

The region of the mass-lifetime plane excluded by this consideration is shown in Fig. 3 
(labeled SNL). 

3.2 Decays outside the progenitor 

The differential gamma-ray fluence at Earth due to tau-neutrino decays via the channel 
v T — > z/ e 7 was calculated in Ref. [^]. As modified to our approximations it is 

dF± = B^K, /rj.y rt- dteM _ 2E t/muTv) ( x + eM _ Emm/Ti/) , (16) 

aE 7 T u m u r u \r rc iJ Jo \ T u J 

where r ns and T v are the radius and temperature of the tau neutrinosphere respectively, r re \ 
is the radius of the light neutrinosphere, E m i n is the maximum of E 1 and 



m,. 



100 sec-m u /2E 7 t) 



1/2 



ciiici i max — 223. 2 sec is the time interval over which the authors of Ref. |7J analysed the data 
recorded by the Gamma-Ray Spectrometer on board the Solar Maximum Mission satellite. 
Eq. ( |i~6t ) differs from the original expression in Ref. in two respects: First, the neutri- 
nosphere temperature T u depends upon tau-neutrino mass; and second, the tau-neutrino 
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fluence is must be modified by the geometrical factor (r ns /r rc i) 2 relative to the light neutrino 
fluence F v = 1.4 x 10 10 cm -2 . 

In this case we are only interested in decays that produce a gamma ray. However, for 
MeV-mass tau neutrinos the mode u T — > u + e ± is kinematically allowed and likely dominates 
the photonic mode. Assuming it proceeds via standard electroweak interactions, 

19271-3 „ lsec 

Tu ~ G 2 F ml sin 2 6 cos 2 9 ~ sin 2 20(m J/ /lO MeV) 5 ' ( ' 

where 9 is the v T — v e mixing angle. To wit, we have taken the factor B 1 ~ 10 -3 , corresponding 
to the fraction of decays expected to have a bremsstrahlung photon (also see Ref. |21|j ). 



Integrating Eq. QIBD over the suitable energy bands and comparing with the corresponding 
3cx fluence limits reported in Ref. leads to the excluded region denoted SMM in Fig. 3 
H2C 



3.3 SN 87A neutrinos 

The nineteen neutrino events recorded by the KII and 1MB detectors confirmed the standard 
model of the cooling of a nascent neutron star and had an energy distribution consistent with 
a temperature of about 4 MeV |Tl|. As noted earlier and illustrated in Figs. 2b and 2c, for 



very short lifetimes (r u < 10~ 6 sec), decays and inverse decays become so potent that they 
move the tau and electron neutrinospheres outward to a region where the temperature is 
around 1 MeV or less. Clearly this is inconsistent with the KII and 1MB data. However, 
we hasten to add that our model for the neutrinosphere is not self consistent under such 
extreme conditions. More importantly, lifetimes for the modes v T — > v e + or v T — » v + 7 
shorter than about 1 sec are excluded by data from the Big European Bubble Chamber 
(BEBC) experiment |19j (shown in Fig. 3). Unless new interactions beyond the standard 
model can give rise to a very short lifetime for a nonelectromagnetic mode where all the 
daughter products would be in thermal equilibrium in a hot neutron star, it would seem that 
our results for very short lifetimes are without application. 

In passing we mention that if tau neutrinos decay beyond the electron neutrinosphere, 
corresponding to r v > 10~ 3 sec, there is an additional population of antielectron neutri- 
nos, which for m u < 10 MeV is of comparable importance. It appears unlikely that our 
understanding of SN 87A and the KII/IMB data are good enough to exclude this possibil- 
ity. However, because MeV-mass neutrinos are only semi-relativistic, these additional events 
would be characterized by a time delay/spread on the order of the tau-neutrino lifetime itself 
and thus could perhaps be excluded (or confirmed) on that basis. (In fact, it could be wildly 
speculated that a tau neutrino with a lifetime of seconds and a mass of 1 MeV — 10 MeV 
could be responsible for the second cluster of KII events.) 

3.4 Residual annihilations 

Finally, we would like to draw attention to an interesting aspect of the behaviour of the 
numerical solutions of Eq. flT0| ) in the long-lifetime limit. In Figs. 4 we show the solution and 
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the equilibrium distribution for the massless case and for m v = 30 MeV. The departure from 
equilibrium near the neutrinosphere radius is clearly visible. The kinks mark the surface of 



last scattering (r = r&s) where the boundary condition Eq. (11) applies. In the massless 
case the actual tau-neutrino number density "undershoots" the equilibrium density in the 
diffusion region, whereas in the massive case it "overshoots." The "overshoot effect" suggests 
that a massive neutrino could be effective in depositing energy into the outer core region by 
residual annihilations and thus could help power the supernova explosion itself. (At present, 



it is not understood precisely how the explosion is powered [11, 17 ). 

To study this further we calculated the energy deposited in the region between the light 
neutrinosphere and the outer boundary of the diffusion region. (Inside the light neutri- 
nosphere any energy deposited will be radiated in electron/muon neutrinos, and beyond the 
diffusion region annihilations are totally ineffective because of the free streaming of neutri- 
nos). We approximate this energy by 

E vD ~ \m v + 3T U ) / Anr dra uP v [n - n eq ) , (18) 

where a u pv is evaluated at the tau neutrinosphere and n(r) is the analytical solution of 
Eq. ([Top between r ns and r^is for vanishing right-hand side and D(r) = D(r ns )(r /r ns ) k . In 
Fig. 5 we show the ratio of the energy deposited to the energy carried off by a massless neu- 
trino species (about 10 53 erg). For tau-neutrino masses between about 15 MeV and 25 MeV 
the energy deposited is about 10 51 erg, which is comparable to the energy seen in the ki- 
netic motion of the explosion and thus could have an important effect on the supernova 
explosion — perhaps solving the riddle of the explosion itself. 



4 Concluding Remarks 

We have studied the transport of energy by MeV-mass tau neutrinos during the early cooling 
phase of a hot neutron star. By means of a simple, but accurate, diffusion model we have 
computed the energy fluence and neutrinosphere temperature. Because the neutrinosphere 
temperature increases with tau-neutrino mass the energy carried off by massive tau neutrinos 
does not decrease nearly as rapidly as it would if the tau-neutrinosphere temperature were 
independent of mass, as previously assumed. (We mention that our method could easily be 
generalized to calculate the flux of hypothetical, massive particles with interactions whose 
strength is roughly weak — e.g., new particles predicted in supersymmetric models — if the 
motivation should arise.) 

Using these new results we have re-derived the mass-lifetime bounds that follow from SN 
87A. Based upon SMM gamma-ray data we exclude a massive tau neutrino of lifetime shorter 
than about 10 8 sec whose decays are radiative. This, taken with the recent bounds based 
upon primordial nucleosynthesis [TR E3], excludes a tau neutrino more massive than about 



0.4 MeV regardless of its lifetime, provided only that its decays are radiative. In addition, 
a Dirac tau neutrino of mass from about lOkeV to 1 MeV is excluded on the basis of the 
fact that wrong-helicity tau-neutrinos produced deep in the core would have carried off too 
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much energy, leading to an unacceptably short burst of antielectron neutrinos . Finally, 
we speculate that a tau neutrino of mass 15 MeV — 25 MeV could play an important role in 
the supernova explosion process itself by virtue of the energy its annihilations deposit just 
outside the tau neutrinosphere. 
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Figure Captions 



Figure 1: Tau-neutrinosphere temperature as a function of tau-neutrino mass: (a) Long- 
lifetime limit (upper curves) and r„ = 10 _5 sec (lower curves) for Dirac (solid) and Majorana 
(broken) mass; (b) Contours in the mass-lifetime plane for Dirac mass; and (c) Contours in 
the mass-lifetime plan for Majorana mass. Contours are in steps of 0.5 MeV beginning with 
7.5 MeV in the upper left corner. The temperature increases to the right and decreases to 
the bottom. 

Figure 2: Massive tau-neutrino energy flux relative to a massless tau neutrino: (a) Long- 
lifetime limit (upper curves) and r u = lCT 5 sec (lower curves) for Dirac (solid) and Majorana 
(broken) mass compared to the naive Boltzmann suppression factor for mass-independent 
neutrinosphere temperature T v = 7.5 MeV (dotted); (b) Contours in the mass-lifetime plane 
for Dirac mass; and (c) Contours in the mass-lifetime plane for Majorana mass. Flux contours 
decrease by a factor of 10 (per contour) from top left to the lower right. 

Figure 3: Excluded regions of the tau-neutrino mass-lifetime plane: (a) Dirac mass; and 
(b) Majorana mass. Forbidden regions are on the same side as the label; SNL is the limit 
based upon decays inside the progenitor; SMM is the limit based upon decays outside the 
progenitor and the excluded region is between the curve just above SMM and the unmarked 
lower curve; CLEO/ARGUS refers to the laboratory mass limit; and BEBC refers to the 
laboratory lifetime limit discussed in Sec. 3.3. 

Figure 4: Numerical solution for tau-neutrino number density (solid curve) and equilibrium 
number density (broken curve), both relative to the neutrino number density at the inner 
boundary (see Sec. 2.3): (a) massless tau neutrino; and (b) 30 MeV Dirac tau neutrino. 

Figure 5: Energy deposited by residual annihilations of a tau neutrino as a function of tau- 
neutrino mass in units of 10 53 erg; solid curve is for Dirac and broken curve is for Majorana. 
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